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Critical properties of the compact three-dimensional U{1) lattice gauge theory are explored at 
finite temperatures. The critical point of the deconfinement phase transition, critical indices and 
the string tension are studied numerically on lattices with temporal extension N =E and spatial 
extension ranging from L ^ 32 to L ~ 256. The critical indices, which govern the behaviour 
across the deconfinement phase transition, are generally expected to coincide with the critical 
indices of the two-dimensional model. It is found that the determination of the infinite volume 
critical point differs from the pseudo-critical coupling at L = 32, found earlier in the literature and 
implicitly assumed as the onset value of the deconfined phase. The critical index v computed from 
the scaling of the pseudocritical couplings agrees well with the value v = 1/2 of the XY model. 
The computation of the index rj brings to a value larger than expected. The possible reasons for 
such behaviour are discussed. 
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1. Introduction 

This article deals with the compact three-dimensional (id) U{\) lattice gauge theory (LGT), 
whose partition function can be written as 

xeA n=Q 

where A is an L? x A'^ lattice, S is the Wilson action, which reads 

S[co] =j8,^cosw(As)+i3f£cosw(A) (1.2) 

ft Pi 

and sums run over all space-like (p^) and time-like (pt) plaquettes. The plaquette angles (o{p) are 
defined in the standard way. The anisotropic couplings pt and jS., are defined in Ref. To study 
the theory at finite temperature, periodic boundary conditions in the temporal direction are imposed 
on the gauge fields. 

At zero temperature the theory is confining at all values of the bare coupling constant 
while at finite temperature the theory undergoes a deconfinement phase transition. It is well known 
that the partition function of the 3d U{1) LGT in the Villain formulation coincides with that of the 
2d XY model in the leading order of the high-temperature expansion [^. When combined with the 
universality conjecture by Svetitsky-Yaffe |Q, this result leads to conclude that the deconfinement 
phase transition belongs to the universality class of the 2d XY model, which is known to have 
Berezinskii-Kosterlitz-Thouless (BKT) phase transition of infinite order [||, ^]. In particular, one 
might expect the critical behaviour of the Polyakov loop correlation function r{R) to be governed 
by the following expressions 

nR) - ^ , (1.3) 

for j8 > and 

r{R) X exp[-R/^{t)] , (1.4) 

for j3 < j3c, t = j8c/j8 — 1. Here, /? » 1 is the distance between test charges, T is the temperature 
and ^ r^e"'-' is the correlation length. Such behaviour of ^ defines the so-called essential scaling. 
The critical indices Tj(r) and v are known from the renormalization-group (RG) analysis of the 
XY model: T] (Tc) = 1/4 and v = 1 /2, where Tc is the BKT critical point. 

The direct numerical check of these predictions was performed on lattices x Nt with L = 
16, 32 and A'^f = 4, 6, 8 in Ref. [0]. Though the authors of Ref. |Q] confirm the expected BKT nature 
of the phase transition, the reported critical index is almost three times that predicted for the XY 
model, T] (Tc) « 0.78. More recent numerical simulations of Ref. [||] have been mostly concentrated 
on the study of the properties of the high-temperature phase. In these papers it was found that, for 
the isotropic lattice jS., = j3, = j8 with L = 32 and A^, = 8, the pseudo-critical point is jSpc = 2.30(2) 
for Ref. [0] and (5^ ~ 2.346(2) for Ref. Values of j3 above these values were taken implicitly 
as belonging to the deconfined phase. 

In Ref. [|I|] we have studied the model on extremely anisotropic lattice with jS^ = 0. By a 
simple analytical analysis we showed that in the limits of both small and large jSj such anisotropic 
model reduces to the 2d XY model with some effective couplings. Then we performed numerical 
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simulations of the effective spin model for the Polyakov loop which can be exactly computed in the 
limit = 0. We used lattices with A^, = 1,4,8 and with the spatial extent L G [64,256] and found 
that the index T] is well compatible with the XY value. We may thus assume that, at least in the 
limit j8,v = 0, the 'idU{\) LGT does belong to the universality class of the XY model. 

Here we consider the isotropic model on the lattice with Nt = 8. Our strategy is the following: 
we postulate that the scaling laws of the XY model are valid and use them to determine the critical 
indices of the gauge model. In doing so we have encountered certain surprises: (i) the infinite 
volume critical coupling turned out to be essentially higher than the values for the pseudo-critical 
couplings reported in Refs. §]. As a consequence, the values of j8 used in Ref. ^ to study 
the deconfinement phase lie well inside the confinement phase when the thermodynamic limit is 
considered; (ii) the index v extracted from the scaling of the pseudo-critical couplings with L 
does agree well with the expected XY value v = 1/2, but the index Tj was found to be strikingly 
different from the XY value, namely r\ 0.50. While the value tj 0.78 obtained in Ref. [Q] 
could, in principle be attributed to rather small lattices used, L = 32, and to an incorrect location of 
the critical point, our result is almost insensitive to varying the spatial extent if L is large enough. 

2. Numerical results 

We simulated the system on lattices of the type x Nt, with Nt = ^ fixed and L increasing 
towards the thermodynamic limit (for details, see Ref. [§]). In the adopted Monte Carlo algorithm 
a sweep consisted in a mixture of one Metropolis update and five microcanonical steps. Measure- 
ments were taken every 10 sweeps in order to reduce the autocorrelation and the typical statistics 
per run was about 100k. The error analysis was performed by the jackknife method over bins at 
different blocking levels. 

The observable used as a probe of the two phases of the finite temperature 3d U{1) LGT is the 
Polyakov loop, defined as 

Pix)=llUoix,t) , (2.1) 

t 

where Uo{x,t) is the temporal link attached at the spatial point x. The effective theory for the 
Polyakov loop is two-dimensional and possesses global U{1) symmetry. Since the global sym- 
metry cannot be broken spontaneously in two dimensions owing to the Mermin-Wagner-Coleman 
theorem, the expectation value of the Polyakov loop vanishes in the thermodynamic limit. On a 
finite lattice (Ly^(-^)) = due to U{1) symmetry (if the boundary conditions used preserve the 
symmetry). This is confirmed by the numerical analysis on the periodic lattice: in the confined 
(small j8) phase the values taken by the Polyakov loop in a typical Monte Carlo ensemble scatter 
around the origin of the complex plane forming a uniform cloud, whereas in the deconfined (high 
j8) phase they distribute on a ring, the thermal average being equal to zero in both cases. What 
really feels the transition is then the absolute value of P, which has been chosen to be the order 
parameter in this work. 

At finite volume the transition manifests through a peak in the magnetic susceptibility of the 
Polyakov loop, defined as 

XL = L\{\P\')-{\P\f) , P=7^lP(^) . (2.2) 
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The value of the couphng at which this happens is the pseudo-critical coupling, jSpc. By increasing 
the spatial volume, the position of the peak moves towards the (nonuniversal) infinite volume crit- 
ical coupling, jSc- The value of jS^c for a given L is determined by interpolating the values of the 
susceptibility Xi around the peak by a Lorentzian function. In Table [I] we summarize the resulting 
values of jSpc and the peak values of the susceptibility for the several volumes considered in this 
work (we included also the determination for L = 32, taken from the first paper in Ref. [||]). 



Table 1: fipc and peak value of the Polyakov loop susceptibility Xi on the lattices x 8. 



L 




^L.max 


32 


2.346(2), Ref. [§ 




48 


2.4238(67) 


12.93(41) 


64 


2.4719(39) 


20.09(66) 


96 


2.5648(96) 


38.8(1.6) 


128 


2.6526(59) 


60.1(3.5) 


150 


2.68(1) 


92.6(8.0) 


200 


2.7336(69) 


144(12) 


256 


2.7780(40) 


220(20) 



There are, in principle, two hypotheses to be tested in order to locate the infinite volume 
critical coupling j3c: first order and BKT transition. The hypothesis of first order transition is not 
incompatible with data for the peak susceptibility for L > 128. However, the corresponding scaling 
law for the pseudo-critical couplings. 



(2.3) 



seems to be ruled out by our data (;^^/d.o.f equal to 5.6 for L > 96, 3.7 for L > 128, 2.1 for L > 96). 

Assuming the essential scaling of the BKT transition, i.e. ^ ^ e'^' ^ 
becomes 

Ppc=Pc + - 



\ the scaling law for jSpc 
(2.4) 



(lnL + B)v 

A 4-parameter fit of the data for Ppc{L) given in Table |l| with the law given in Eq. ( ^4| ) leads to 
unstable values of the parameters. Instead, when the parameter v is fixed at the XY value, v = 1/2, 
the fit is stable for lattices with size not smaller than L, leading to an estimated value of the infinite 
volume critical couphng, jSc = 3.06(11) (see Ref. |^). 

By finite size scaling (FSS) analysis at jSc, we can extract other critical indices. An interesting 
one is the magnetic critical index, Tj , which enters the scaling law 



'2-11 



(2.5) 



Actually in this law one should consider logarithmic corrections (see Refs. [|10[ [11|] and references 
therein) and, indeed, recent works on the XY universality class generally include them. However, 
taking these corrections into account for extracting critical indices calls for very large lattices even 
in the XY model; for the theory under consideration to be computationally tractable, we have no 
choice but to neglect logarithmic corrections. 
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Setting the coupling j3 at the value of our best estimation for j3c, i.e. j3 = 3.06, we determined 
the susceptibilities XiiPc) for several volumes. Then, following FSS, we fitted the results with the 
law XiiPc) =AL^^'^ and got 

A = 0.0171(10), T7 = 0.496(15) (x Vd-o.f. = 0.60) . (2.6) 

This value for rj is by far incompatible with the 2d XY value, rjxY = 0.25. The most ex- 
treme consequence of this finding is that the deconfinement transition in the 3d U{1) LGT at finite 
temperature does not belong to the same universality class as 2d XY spin model. This would con- 
tradict the Svetitsky-Yaffe conjecture, raising a problem in the understanding of the deconfinement 
mechanism in gauge theories. We will further comment on this issue in the discussion section. 

In such a situation, it becomes particularly useful to have another determination of the index 
77, by an independent approach. Following Ref. we define an effective rj index, through the 
2-point correlator of Polyakov loops, according to 

„ , . _ iog[r(j?)/r(/?o)] 

log[Ro/R\ 

with/?o chosen equal to 10, as in Ref. This quantity is constructed in such a way that it exhibits 
a plateau in R if the correlator obeys the law ([T^, valid in the deconfined phase. 

The analysis of the behaviour of T]eff(/?) has been repeated setting j3 at our estimated value 
for j8c, i.e. j3 = 3.06, and increasing the spatial extent of the lattice. It turns out (see Fig. [I]) that 
a plateau develops at small distances when L increases and that the extension of this plateau gets 
larger with L, consistently with the fact that finite volume effects are becoming less important. 
The plateau value of T]eff can be estimated as T]eff(/? = 6) on the 256^ x 8 lattice and is equal to 
0.4782(25); it agrees with our previous determination of the index r]. 



3. Discussions 

We have studied the critical behaviour of the 3dU{\) LGT at finite temperature on isotropic 
lattice with the temporal extension Nt = S. The pseudo-critical coupling was determined through 
the peak in the susceptibility of the Polyakov loop; the infinite-volume critical coupling has then 



been computed assuming the scaling behaviour of the form (p^), the result being jSc = 3.06(11). 
The deconfinement phase is the phase where j8 > jS^. A thorough investigation of the deconfine- 
ment phase was performed in Ref. [||]. However, all j8 -values used there are smaller than the 
infinite-volume critical coupling. When the thermodynamic limit is approached the critical cou- 
pling increases so that the numerical results of Ref. ^ would refer rather to the confinement phase 
of the infinite-volume theory. 

We found also that the index 17 turns out to be T] « 0.496. This value is essentially larger than 
expected and requires some discussion. The easiest explanation would be to state that the spatial 
lattice size used (L G [32 — 256]) is still too small to exhibit the correct scaling behaviour, hence 
the wrong values for jSc and rj follow. However, if one makes a plot of Ppc{L) vs L, one can see, by 
looking at the trend of data, that it is unlikely that j8c is much larger than our estimate. In fact, our 



fits with the scaling law (2.4) show that pc decreases when larger lattices are considered. Therefore, 
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Figure 1: T7ei^ vs R at j5c on lattices with several values of L. 



our result is most likely an overestimation. This implies that the true 7] is likely even larger than 
what we found. 



Moreover, neglecting logarithmic corrections to the scaling law ( |2.5[ ) cannot have such a strong 
impact to decrease 77 to half of the value we found. 

Let us give a simple argument why the index 77 can be different from its XY value. Consider the 
anisotropic lattice, in the limit of large jS^ = °°. Here the spatial plaquettes are frozen to unity and the 
ground state is a state where all spatial fields are pure gauge, i.e. Un{x) = VtV/+^„ , n = 1,2. Under 
the change of variables Uo{x) — > VxUo{x)V*ji_gg, in the leading order of the large- jS^ expansion the 
partition function factorizes into the product of A^, independent 2d XY models. Since the Polyakov 
loop is the product of gauge fields in the temporal direction, the correlation function factorizes, too, 
and becomes a product of independent XY correlations, i.e. 



r[/(i)(Av = -,i8r) 

Hence, for asymptotically large /? ^ 1 , we get 

rt/(i)(i8. = oo,ft>r) 



RVxy 



N, 



(3.1) 



(3.2) 



This leads to a simple relation 

T](As.=oo,j3r) = A^rT]xF. (3.3) 

Some conclusions could now be drawn. The critical behaviour of the 3dU{l) LGT in the limit 
^ 00 is also governed by the 2d XY model. Nevertheless, the effective index 77 appears to be Nt 



6 



Critical behaviour of the compact 3dU{l) gauge theory at finite temperature 



Roberto Fiore 



times of its XY value. Now, for jS^ = we have r\[^s = ^,PD = Vxy- This relation and formula 



( |3.3| ) allow to conjecture that 

rixY < Tli^,Pn < NtilxY ■ (3.4) 

Ps = corresponds to the lower limit while j8s = oo corresponds to the upper limit. In general, 
77 could interpolate between two limits with p^- Whether this interpolation is monotonic or there 
exists critical value jS/'', such that 77(j8.v < j8/'',j3"") = rjxy and 77 changes monotonically above 
P", cannot be answered with data we have and requires simulations on the anisotropic lattices. In 



the paper [12| a renormalization group study of 3d U{\) model at small jS^ will be presented and 
computations of the leading correction to the large jS^ behaviour will be given. The results of our 
computations support the scenario that the index rj depends on the ratio Ps/Pt- Recently, we have 



obtained the results of simulations for Nt = 2,4 performed by A. Bazavov []13|]. His results also 
point in the direction of our scenario. 

Finally, it is worth mentioning that the factorization in the large jS., limit does not affect the 
index v. It follows from its definition ( |1.4| ) that in this limit v = 1/2 as in the XY model. We expect 
therefore that v equals 1 /2 for all p^ and is thus universal. 

In view of our results it might be worth to perform numerical simulations for small but nonva- 
nishing ps and for larger volumes. 
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